On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

A high precision U–Pb age of metamorphic rutile in coesite-bearing eclogite from the Dabie Mountains in central China: a new constraint on the cooling history

This paper first reports a high precision U–Pb age of 218±1.2 Ma for rutile in coesite-bearing eclogite from Jinheqiao in the Dabie Mounteins, east–central China. This work shows that the U–Pb mineral (rutile+omphacite) isochron age of 218±2.5 Ma and conventional rutile U–Pb concordia age of 218±1.2 Ma obtained by common Pb correction based on the Pb isotopic composition of omphacite in the same eclogite sample are consistent, proving that the omphacite with low U/Pb ratio (μ=2.8) can be used for common Pb correction in U–Pb dating of rutile. Oxygen isotope analysis of rutile aliquots gave the consistent δ¹⁸O values of −6.1±0.1%, demonstrating oxygen isotope homogenization in the rutile of different grains as inclusion in garnet and grain in matrix. Oxygen isotope thermometry yields temperatures of 695±35 and 460±15 °C for quartz–garnet and quartz–rutile pairs, respectively. These oxygen isotopic observations suggest that the diffusion of oxygen in rutile as inclusion in garnet is not controlled by garnet. According to field-based thermochronological studies of rutile, an estimate of the T_c of about 460 °C for U–Pb system in rutile under rapid cooling conditions (20 °C/Ma) was advised. Based on this U–Pb age as well as the reported chronological data with their corresponding metamorphic and/or closure temperature, an improved T–t path has been constructed. The T–t path confirms that the UHPM rocks in South Dabie experienced a rapid cooling following the peak metamorphism before 220 Ma and a long isothermal stage from 213 to 180 Ma around 425 °C.     

Lower crustal melting and the role of open-system processes in the genesis of syn-orogenic quartz diorite–granite–leucogranite associations: constraints from Sr–Nd–O isotopes from the Bandombaai Complex, Namibia

The Bandombaai Complex (southern Kaoko Belt, Namibia) consists of three main intrusive rock types including metaluminous hornblende- and sphene-bearing quartz diorites, allanite-bearing granodiorites and granites, and peraluminous garnet- and muscovite-bearing leucogranites. Intrusion of the quartz diorites is constrained by a U–Pb zircon age of 540±3 Ma.

Quartz diorites, granodiorites and granites display heterogeneous initial Nd- and O isotope compositions (_{Nd (540 Ma)}=−6.3 to −19.8; δ¹⁸O=9.0–11.6‰) but rather low and uniform initial Sr isotope compositions (⁸⁷Sr/⁸⁶Sr_initial=0.70794–0.70982). Two leucogranites and one aplite have higher initial ⁸⁷Sr/⁸⁶Sr ratios (0.70828–0.71559), but similar initial _Nd (−11.9 to −15.8) and oxygen isotope values (10.5–12.9‰). The geochemical and isotopic characteristics of the Bandombaai Complex are distinct from other granitoids of the Kaoko Belt and the Central Zone of the Damara orogen. Our study suggests that the quartz diorites of the Bandombaai Complex are generated by melting of heterogeneous mafic lower crust. Based on a comparison with results from amphibolite-dehydration melting experiments, a lower crustal garnet- and amphibole-bearing metabasalt, probably enriched in K₂O, is a likely source rock for the quartz diorites. The granodiorites/granites show low Rb/Sr (<0.6) ratios and are probably generated by partial melting of meta-igneous (intermediate) lower crustal sources by amphibole-dehydration melting. Most of the leucogranites display higher Rb/Sr ratios (>1) and are most likely generated by biotite-dehydration melting of heterogeneous felsic lower crust. All segments of the lower crust underwent partial melting during the Pan-African orogeny at a time (540 Ma) when the middle crust of the central Damara orogen also underwent high T, medium P regional metamorphism and melting. Geochemical and isotope data from the Bandombaai Complex suggest that the Pan-African orogeny in this part of the orogen was not a major crust-forming episode. Instead, even the most primitive rock types of the region, the quartz diorites, represent recycled lower crustal material.     

Accretion and exhumation at a Variscan active margin, recorded in the Saxothuringian flysch

The Saxothuringian flysch basin, on the north flank of the Central European Variscides, was fed and eventually overthrust by the northwestern, active margin of the Tepla-Barrandian terrane. Clast spectra, mineral composition and isotopic ages of detrital mica and zircon have been analyzed in order to constrain accretion and exhumation of rocks in the orogenic wedge. The earliest clastic sediments preserved are of early Famennian age (ca. 370?Ma). They are exposed immediately to the NW of the suture, and belong to the par-autochthon of the foreland. Besides ultramafic (?ophiolite) material, these rocks contain clasts derived from Early Paleozoic continental slope sediments, originally deposited at the NW margin of the Saxothuringian basin. These findings, together with the paleogeographic position of the Famennian clastics debris on the northwestern passive margin, indicate that the Saxothuringian narrow ocean had been closed by that time. Microprobe analyses of detrital hornblendes suggest derivation from the “Randamphibolit” unit, now present in the middle part of the Saxothuringian allochthon (Münchberg nappes). Detrital zircons of metamorphic rocks formed a little earlier (ca. 380?Ma) indicate rapid recycling at the tectonic front. The middle part of the flysch sequence (ca. early to middle Viséan), both in the par-autochthon and in the allochthon, contains abundant clasts of Paleozoic rocks derived from the northwestern slope and rise, together with debris of Cadomian basement, 500-Ma granitoids and 380?Ma (early Variscan) crystalline rocks. All of these source rocks were still available in the youngest part of the flysch (c. middle to late Viséan), but some clasts record, in addition, accretion of the northwestern shelf. Our findings permit deduction of minimum rates of tectonic shortening well in excess of 10–30?mm per year, and rates of exhumation of ca. 3?mm/a, and possibly more.     

Sampling Fauna in Stream Sediments as well as Groundwater Using One Net Sampler

Groundwater bores act as traps. Net samplers are regularly used for sampling this type of trap for fauna. To enable direct comparisons of faunal communities in groundwater bores and stream sediments, stream sediment tubes were built similar to groundwater bores and were sampled with net samplers for fauna. These stream sediment tubes consisted of a tube anchored in the stream sediment, also called interstitial space. To test the efficacy of this trap method in stream sediments, it was compared to another type of trap, Hahn's trap. Faunal communities sampled by a net in the stream sediment tubes did not differ hugely from fauna in Hahn's trap samples. Physical and chemical factors of sampled water in both the stream sediment tubes, the surrounding interstitial sediments and the second type of traps, Hahn's traps, showed that water in both the tubes and Hahn's traps was closely related to interstitial water. The net sampler is inexpensive and easy to handle. It is suggested that sampling stream tubes with nets may be an appropriate method for long‐term monitoring studies.     

Modelling transverse dunes

Transverse dunes appear in regions of mainly unidirectional wind and high sand availability. A dune model is extended to two‐dimensional calculation of the shear stress. It is applied to simulate dynamics and morphology of three‐dimensional transverse dunes. In the simulations they seem to reach translational invariance and do not stop growing. Hence, simulations of two‐dimensional dune ?elds have been performed. Characteristic laws were found for the time evolution of transverse dunes. Bagnold's law of the dune velocity is modi?ed and reproduced. The interaction between transverse dunes led to the interesting conclusion that small dunes can travel over bigger ones. Copyright © 2004 John Wiley & Sons, Ltd.     

Kinetic energy–rainfall intensity relationship for Central Cebu, Philippines for soil erosion studies

In the study of soil erosion, specifically on detachment of soil particles by raindrop impact, kinetic energy is a commonly suggested indicator of the raindrop's ability to detach soil particles from the soil mass. Since direct measurement of kinetic energy requires sophisticated and costly instruments, the alternative approach is to estimate it from rainfall intensity. The present study aims at establishing a relationship between rainfall intensity and kinetic energy for rainfalls in Central Cebu, Philippines as a preface of a wider regional investigation.

Drop size distributions of rainfalls were measured using the disdrometer RD-80. There are two forms of kinetic energy considered here. One is kinetic energy per unit area per unit time (KE_R, J m⁻² h⁻¹) and the other is kinetic energy per unit area per unit depth (KE, J m⁻² mm⁻¹). Relationships between kinetic energy per unit area per unit time (KE_R) and rainfall intensity (I) were obtained using linear and power relations. The exponential model and the logarithmic model were fitted to the KE–I data to obtain corresponding relationships between kinetic energy per unit area per unit depth of rainfall (KE) and rainfall intensity (I). The equation obtained from the exponential model produced smaller standard error of estimates than the logarithmic model.     

Denitrification rates measured along a salinity gradient in the eutrophic Neuse River estuary, North Carolina, USA

Denitrification rates along a salinity gradient in the eutrophic Neuse River Estuary, North Carolina, were quantified using membrane inlet mass spectrometry (MIMS) within short-term batch incubations. Denitrification rates within the system were highly variable, ranging from 0 to 275 μmol N m⁻² h⁻¹. Intrasite variability increased with salinity, but no significant differences were observed across the salinity gradient. Denitrification rates were positively correlated with sediment oxygen demand at the upstream sampling site where sediment organic carbon levels were lowest. This relationship was not observed in the more saline sampling sites. Denitrification rates were highest during winter. On an annual basis, denitrification accounted for 26% of the dissolved inorganic nitrogen and 12% of the total nitrogen supplied to the system.     

Importance of climate,forest fires and human population size in the Holocene boreal forest composition change in northern Europe

The relative importance of climate, forest fires and human population size on long‐term boreal forest composition were statistically investigated at regional and local scales in Fennoscandia. We employ pollen data from lakes, reflecting regional vegetation, and small forest hollows, reflecting local vegetation, from Russia, Finland and Sweden to reconstruct the long‐term forest composition. As potential drivers of the Holocene forest dynamics we consider climate, generated from a climate model and oxygen isotope data, past forest fires generated from sedimentary charcoal data and human population size derived from radiocarbon dated archaeological findings. We apply the statistical method of variation partitioning to assess the relative importance of these environmental variables on long‐term boreal forest composition. The results show that climate is the main driver of the changes in Holocene boreal forest composition at the regional scale. However, at the local scale the role of climate is relatively small. In general, the importance of forest fires is low both at regional and local scales. The fact that both climate and forest fires explain relatively small proportions of variation in long‐term boreal vegetation in small forest hollow records demonstrates the complexity of factors affecting stand‐scale forest dynamics. The relative importance of human population size was low in both the prehistorical and the historical time periods. However, this is the first time that this type of data has been used to statistically assess the importance of human population size on boreal vegetation and the spatial representativeness of the data may cause bias to the analysis.